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Computational  Fluid  Dynamics  (CFD):  Element-by-Element 

Analysis  ^ 

for  Advanced  Computers  THirV 

Principal  Investigator:  G.  F.Carey,  Director  -  v  ' 
CFD  Laboratory,  University  of  Texas 


Overview;  The  research  has  focussed  on  a  class  of  domain  decomposition 
techniques  for  solution  of  problems  in  Computational  Fluid  Dynamics  using 
advanced  vector  and  parallel  processors.  In  particular  we  consider  finite 
element  schemes  and  use  the  natural  element-partition  of  the  domain  to 
constmct  the  decomposition  algorithm.  This  then  fits  conveniently  into  the 
usual  framework  of  finite  element  calculations  in  which  the  primary  loop  is 
the  independent  calculation  and  assembly  of  element  matrix  and  vector 
contributions.  By  recasting  the  conjugate  gradient  method  at  this  level,  the 
system  matrix  need  not  be  assembled  and  the  intensive  matrix-vector  product 
step  can  be  completely  parallelized.  These  ideas  constitute  a  major  departure 
from  traditional  finite  element  schemes  and  we  feel  our  efforts  are  a  major 
new  development  that  will  strongly  influence  the  technology.  One  difficulty 
for  CFD  applications  is,  however,  the  convective  term  that  leads  to  highly 
non-symmetric  systems.  The  conjugate  gradient  scheme  will  not  converge 
for  such  systems.  To  circumvent  this  problem  we  have  been  investigating 
and  implementing  biconjugate  gradient  and  conjugate  gradient  squared 
techniques.  Preliminary  results  indicate  that  we  have  resolved  this  difficulty. 
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Status:  We  have  analyzed  and  implemented  the  EBE  scheme  on  the  CRAY 
XMP  at  Texas,*  and  developed  a  vectorized  scheme  with  vectorization  over 
the  number  of  elements.  Hence  vector  lengths  are  very  long.  This  scheme  is 
also  being  implemented  on  the  ETA- 10  at  Florida  but  we  are  still  "ironing 
out"  software  problems.  The  algorithm  for  parallel  processing  CFD  is 
different  and  has  been  examined,  thus  far,  for  parallel  scalar  shared- 
memory  machines.  Our  calculations  of  test  problems  have  been  made  at 
Argorme  on  the  Sequent  Balance  (20  processors  in  parallel)  and  the  Alliant 
FX-8  (8  processors  in  parallel).  In  late  July  we  acquired  an  Alliant  FX-8  at 
TICOM  and  the  CFD  lab  is  connected  via  ethemet  to  the  Alliant  (  1  floor 
apart  on  the  Dept,  ethemet  loop).  We  are  now  using  this  local  Alliant  and 
this  will  expedite  our  research  and  application  studies. 

Related  Activities 

ADVANCED  PROCESSOR  SYMPOSIUM:  On  Oct.  17, 18 1  organized  and 
hosted  a  symposium  on  "Methods  and  Algorithms  for  PDE’s  on  Advanced 
Processors".  A  copy  of  the  program  and  list  of  speakers  is  attached. 

ITERATIVE  METHODS  CONFERENCE:  A  conference  on  Iterative 
Methods  was  organized  with  support  from  DOE,  AFOSR,  and  with  the 
sponsorship  of  SIAM,  CNA  and  the  CFD  Lab.  Over  150  attendees  from  the 
US  and  overseas  participated  in  the  3-day  conference  which  was  organized  to 
honor  Dr.  Young's  65th  birthday  and  longstanding  contributions  to  this  field. 
All  aspects  of  the  conference  were  a  great  success.  A  copy  of  the  conference 
program  is  attached.  A  conference  volume  of  invited  papers  is  to  be 
prepared  and  contributed  papers  will  appear  in  the  journals  Communications 


in  Applied  Numerical  Methods  and  International  Journal  of  Numerical 
Methods  in  Engineering. 
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Reprints  and  Preprints 

Element-by-Element  Vector  and  Parallel  Computations,  CANM,  4,  No.  3, 
299-307, 1988  (with  E.  Barragy,  R.  McLay  and  M.  Sharma). 

Supercomputing  in  Computational  Fluid  Dynamics,  Computational  Fluid 
Dynamics.  G.  de  Vahl  Davis  and  C.  Fletcher  (Eds.),  Elsevier  Science 
Publishers  B.V.  (North  Holland),  21-30, 1988. 

Least-Squares  Finite  Elements  for  First-Order  Hyperbolic  Systems,  UNME, 
26,  81-93, 1988,  (with  B.  N.  Jiang). 

A  Parallel  Element  by  Element  Solution  Scheme,  UNME  (in  press),  March 
1988,  and  CNA  Report  214  (with  E.  Barragy). 

Extension  of  Inverse  Design  Techniques  for  Multicomponent  Airfoils,  AIAA 
Journal,  (in  press),  1988,  (with  M.  Siladic). 


Presentations 

Keynote  Lecture,  Supercomputing  in  Computational  Fluid  Dynamics, 
Computational  Fluid  Dynamics  Meeting,  University  of  Sydney,  Australia, 
August  1987. 

Distinguished  Massman  Lecture,  Advances  in  Finite  Element  Methods  and 
Supercomputing,  The  University  of  Notre  Dame,  Indiana,  May  1988. 

Seminar,  Recent  Advanced  in  Supercomputing  and  Computational  Fluid 
Dynamics,  Fluids  Group,  The  University  of  Texas  at  Austin,  April  1988. 

In  August,  1987,  Dr.  Carey  delivered  the  keynote  lecture  opening  the  CTAC 
Computational  Techniques  and  Applications  conference  in  Sydney, 
Australia.  The  lecture  will  appear  in  the  bound  conference  proceedings. 


Dr.  Carey  gave  an  invited  lecture  at  Griffiths  University,  Brisbane 
Australia,  August  9  on  "Supercomputers  in  Finite  Element  Computations". 

In  Febmary  1988,  Dr.  Carey  taught  a  short  course  on  computational  aspects 
of  finite  elements  that  included  adaptive  grid  strategies,  element-by-element 
techniques  and  advanced  processors. 


Students 

M.  Sharma,  PhD  Program 
J.  Batista,  MS  Program 
J.  Geld,  MS  Program 
E.  Barragy,  MS  Program 


New  Developments 

Development  and  refinement  of  the  element-by-element  domain 
decomposition  technique  for  Computational  Fluid  Dynamics  on  vector  and 
parallel  computers. 
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